Facilitated spin models in one dimension: a real-space renormalization group study 
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We use a real-space renormalization group (RSRG) to study the low temperature dynamics of ki- 
netically constrained Ising chains (KCICs) . We consider the cases of the Fredrickson- Andersen (FA) 
model, the East model, and the partially asymmetric KCIC. We show that the RSRG allows one 
to obtain in a unified manner the dynamical properties of these models near their zero-temperature 
critical points. These properties include the dynamic exponent, the growth of dynamical length- 
scales, and the behaviour of the excitation density near criticality. For the partially asymmetric 
chain the RG predicts a crossover, on sufficiently large length and time scales, from East-like to 
FA-like behaviour. Our results agree with the known results for KCICs obtained by other methods. 
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I. INTRODUCTION 



II. KCIC: A PSEUDOSPIN FORMULATION 
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Kinetically constrained models (KCMs) [HQ- ES LI are 
systems in which certain trajectories between configura- 
tions are suppressed y|. As a result, they display inter- 
esting slow dynamical behaviour 0, Q, H, 0, E3 ■ Simple 
KCMs, like the facilitated kinetic Ising model introduced 
by Fredrickson and Andersen [| (hereafter the FA model) 
and Jackie and Eisinger Q (hereafter the East model), 
display the slow, cooperative relaxation characteristic of 
supercooled liquids near the glass transition . For 
a general review of KCMs see • 

In this paper we show that a simple real-space renor- 
malization group (RSRG) scheme yields the dynam- 
ical properties of facilitated spin models in one dimen- 
sion, or kinetically constrained Ising chains (KCICs), 
near their zero temperature critical point [l^. This be- 
haviour is known from previous work |l3j . The RSRG 
scheme provides a unified framework for treating these 
systems, allowing one to obtain critical dynamic expo- 
nents, and to visualise the RG flows of scaling variables 
related to temperature and spatial asymmetry. 

We proceed as follows. In Section[H]we define KCICs, 
and show, following Siggia (lS|. that they can be written 
in terms of "interacting quantum spins" . We discuss in 
Section lllll how this formalism admits a simple physical 
interpretation in terms of reaction-diffusion processes. In 
Sections IIVHVIII we use a RSRG scheme to extract the 
zero-temperature critical behaviour of KCICs, in the FA 
model (Sec. EJ and East model (Sec. IVIfl limits, and for 
the case of finite asymmetry |2fJ (Sec lVIljl . We find that 
the dynamical exponents for the FA and East models 
are respectively z = 2 and z = (Tin 2) , in agreement 
with existing results 01 ■ We show that length scales in 
the FA model grow as £ ~ e 1 ^ near the critical point 
T = 0, while for the East model there is no characteristic 
length scale. We also quantify the crossover of the KCIC 
with large but finite asymmetry from East-like to FA-like 
behaviour. In Section IVIIII we state our conclusions. 



The KCIC is defined as follows [HQ. Consider a chain 
of N Ising spins cr, = ±1, in one space dimension, with 
Hamiltonian H = i Yli a^. We will take N even, assume 
periodic boundary conditions, and restrict the dynam- 
ics to flips of single spins that have at least one nearest 
neighbour in the up state. The transition rates depend 
on whether the facilitating up-spin f lies to the left or 
the right of the flipping spin: 



TT^IT, IT^H, H^U, U-^TT, 



(i) 



where b € [0, 1], b = 1 - b, and c = (1 + e}l T y x w e^ 1 ^ 
at low temperature. The bias b determines the sym- 
metry properties of the kinetic constraint: the FA and 
East models correspond to the limiting cases of symme- 
try (& = i) and maximal asymmetry (b — 0), respec- 
tively. The East model is so-called because information 
propagates to the east. We will also consider the case of 
general b, which we will call the biased constrained Ising 
chain (BCIC). For b small but finite the BCIC exhibits a 
crossover at large le ngth and time scales from East-like 
to FA-like behaviour |2(|. In section lTvl we use an RSRG 
to quantify this crossover. 

The dynamics of the KCIC is governed by the master 
equation 



dP{a,t) 
di 



= -J2 w{°i)P{o, t) +]T w(-ai)P(a', t), (2) 



where P(cr, t) is the probability that the system has 
configuration a = {<7i, . . . , <7,, . . . , crjy} at time t (a' is 
the configuration a with spin a% flipped), and w{pi) = 
w((Ti, {<Jj}) is the probability per unit time that Ui will 
flip. The {<7j} are the nearest neighbours of i. The 
matrix controlling the time development of the 2^- 
component vector P is not in general Hermitian, but 
can be made so by introducing the vector ip(a,t) = 

-Po(c) 5 P(v, t). Here Pq(<j) is the equilibrium distribu- 
tion. However, this obscures the fact that the evolution 



2 



operator is a normalized stochastic process which obeys 
detailed balance, and so we will use the non-Hermitian 
representation where this is explicit. 

One passes to a quantum formalism |l9j| by introducing 
the state vector 



where the Liouvillian £, 



1 is 



\P(t)} 



P(a z ,t) |of> <g) |of> ■ 



K> • (3) 



The ket |af) is the state vector for the spin at site i. 
The spins are for convenience taken to lie along the z- 
axis, and thus the operator of flips the z— component 
of spin i; <rff(cr^) = /(— <r? )af. By differentiating © 
with respect to time, and using to eliminate P(a,t), 
we get the master equation in the guise of a Euclidean 
Schrodinger equation, 

^\P(t)) = -H\P(t)). (4) 

Here 7i, which we will also call a "Hamiltonian" , is not 
in general Hermitian. 

We define the KCIC by the constrained Glauber rates 



w(<7i,-0j}) =C l ({crj})- 



o/3cri/2 



(5) 



JJ/ 2cosh(/3/2) 

where the factor of 2 cosh(/3/2) is a convenient normaliza- 
tion. The constraint is Ci{{<jj}) — brii-i + brii+i, where 
rii - 2 

j8 



i(l + cr?). The sigma matrices at a given site i obey 

cfcrj = (5 a/3 + ie"P(Tj. Sigma matrices at different sites 
commute. The matrix 7i in Equation {2J then reads 



(6) 



where TV -1 = 2cosh(/3/2). The appearance of the af 
term shows that in order to represent these simple KCMs 
in terms of "interacting" systems with no dynamical con- 
straints, one must introduce extra degrees of freedom, or, 
equivalently, nonlocal interactions. Using the Pauli rep- 



resentation <r x = 



and a z = 



and 



L = 



1 

1 

5(1 + *f) 





-1 



1 




we have 



(7) 



1 — c — c 
c — 1 c 



(8) 



2cosh/3/2 

The Hamiltonian Tl can then be written as the matrix 
direct product 



JV-l 



H 



(1-6) ^ l®---®ni-i®li 



N-l 

i=l 

N-l 



•iii ® 1 



4=1 



(9) 



£ = {l-b)n®l + M®n 

( 1 - c -fee -be \ 

b(c-l) 6c 

b(c -1)0 be 

\ 0/ 



(10) 



The matrix l|10|l describes a probability-conserving 
stochastic process; thus the sum of each column is zero. 
When we construct the RG scheme we must preserve this 
condition. The East and FA models correspond to the 
cases 6 = and 6 = i, respectively: 



n <8> £ = 



1 - c -c 
c-1 c 





(11) 



and 



1 



-FA 



{n®£ + £(gin) 



(2 -2c 
c-1 
c-1 
V 



In the next section we show briefly that the evolution 
operators obtained above have a simple physical inter- 
pretation as reaction-diffusion processes. 



III. INTERPRETATION VIA 
REACTION-DIFFUSION PROCESSES 

The interpretation of the KCIC as a reaction-diffusion 
process follows by noting that the Liouville operators in 
the previous section act on two-site basis states 



I P(Q) 



P(M) 
P(1,0) 
P(0,1) 
P(0,0) 



(13) 



where P(l) is the probability that a spin is up. We have 
again suppressed time labels. The single-site basis states 

— — 

where p < 1 . The vectors ( J ~\ and ( ^ ^ are the eigen- 
vectors of <7 Z with eigenvalue +1 and —1, respectively. 

We represent reaction-diffusion processes as follows. 
Let the eigenvalue of = 1 represent a lattice site i oc- 
cupied by a particle A, and af = — 1 represent the same 
lattice site with no particle present. We denote this state 
by 0. Then the East model Liouvillian (|llfl has a clear 
physical interpretation: it describes the (right) branch- 
ing process A + 9 — > A + A occurring with rate c, and the 
(right) coagulation process -A + .A — ► A + occurring with 
rate 1 — c. The FA model involves in addition the (left) 
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branching and coagulation processes A + A — > + A and 
Q + A^A + A. 

Consider the following general set of reaction-diffusion 
processes: 



Process 


Description 


Rate 


right diffusion 


A- 


f 0- 


-04 


- A 


Dr 


left diffusion 


0H 


- A - 


■+ A- 


f 


D L 


right coagulation 


A - 


- A- 


-> A 


f 


ICR 


left coagulation 


A - 


- A- 




-A 


1CL 


pair annihilation 


A- 


f A - 


^0- 


f 


1A 


right death 


0- 


- A- 


^0H 


-0 


s 


left death 


A 


+ 0- 


^0H 


-0 


6 


right branching 


A - 


-0- 


* A H 


-A 


PR 


left branching 


+ 


A - 


♦ AH 


-A 


PL 


pair creation 


0H 


-0- 


> AH- 


-A 


V 


right birth 


0- 


r0- 


^0 + 




a 


left birth 


0- 


r0- 




-0 


a 



By inspection, using (|f 3|) , the equation of motion for this 
system reads 



Pit) 



( r -pl -*> \ 

"7ci?, -Di?. —Dl —o 

~1CL ~Dr D L -a 

\ -1A -5 -5 S J 



P(t), (14) 



where T ee 7 Cfl + 7 CL + j A , D R = D R + p R + 5, D L = 
D L + PL + 5 and Y, = 2a + v. 

As an aside, we argued in Rcf. 15] that a d > 1 gen- 
eralization of the FA model behaves like a system in the 
directed percolation (DP) |2l| universality class, albeit 
with vanishing self-destruction rate. The latter process 
may be defined as a reaction-diffusion system compris- 
ing diffusion, branching and annihilation. The similar- 
ity between the FA model and DP is due to the fact 
that nearest neighbour-facilitated branching can mimic 
the effect of diffusion in higher dimensions, and thus the 
d > 1 FA model can be represented as (pseudo)-diffusion, 
branching and annihilation. In one dimension this cor- 
respondence no longer holds: one can see from Equation 
(|14fl that left and right branching, pr and pl, and left 
and right diffusion, Dl and Dr, sit in different slots of 
the evolution matrix. 



is a numerical scheme which tends to produce more pre- 
cise estimates for critical exponents than does the RSRG, 
but is less intuitive, in the sense that it does not lend it- 
self so readily to the visualisation of flows in RG space. 
Field theoretic RG techniques have been successfully ap- 
plied to many reaction-diffusion systems in low dimen- 
sions 0,0. However, we argued in Ref. 0] that the 
FA model has an upper critical dimension of 4, and, cru- 
cially, its coarse-grained action becomes difficult to anal- 
yse below d = 2. Hence we shall employ a real-space RG 
scheme in d = 1, which provides both an intuitive and a 
tractable means of studying the critical behaviour of the 
KCIC. 

The idea is as follows. One divides the lattice into 
blocks of p spins, and denotes the configuration of spins 
inside each block a as {ai} a . We will focus on the case 
p = 2, and discuss larger blockings where appropriate. 
The evolution operator Tip then splits into an intra-block 
part Tio and an inter-block 'interaction' V: 

n = J2( n 0;a+V a , a+1 ). (15) 
a 

In the case of the East model we can write 

Ho = {n <g> £) <g> (1 ® 1), 



and 



V = (1 On) (g) (£ (g) 1), 



(16) 



(17) 



where the brackets indicate the grouping of spins into 
blocks. All terms in the Liouvillian @ are of the form 
(Cnj or (UZJ), with the necessary number of identity ma- 
trices affixed at each end of the chain. 

We will denote the eigenstates of Ho- a as |n) and 
(n|, noting that the left and right eigenstates of a non- 
Hermitian matrix are in general different. The East 
and FA models have respectively triply- and doubly- 
degenerate ground states, i.e. respectively three and two 
eigenvectors with eigenvalue zero. One identifies these 
ground states as effective cell states, and projects the 
basis of two-site spins \a) G {|TT) , |TI) , lit) > III)} mto 
renormalized block spins \a') G {|+) , |— )}. This is done 
by defining a projection operator 



TUaV)^, 



(18) 



IV. REAL-SPACE RENORMALIZATION 
GROUP IN 1 + 1 DIMENSIONS 

We will now apply a simple real-space RG scheme to 
the KCIC. This scheme was developed in the 1980s and 
used on quantum spin models 123. Recently, it was 
applied to the contact process [lij. a reaction-diffusion 
system. Other RG approaches that have been used 
on reaction-diffusion systems include density matrix RG 
(DMRG) [3 and field-theoretic RG techniques. DMRG 



where 



is a linear combination of the the renormalized 



basis vectors \cr'), and the c n , n ' are real numbers. If the 
original Hamiltonian Ti is 2 N dimensional, then T\ is a 
2-/V/2 x 2^ m atrix whose rows contain linear combinations 
of the left ground states of 7io . One can therefore write 
the matrix representation of T\ as 



®N/2 



Ti = 



c ITT),l+> c ITl):l+> c liT).|+> c IU),l+> 
c ITT},|-> c IU},|-> c llT),h> c IU>,|-> 



(19) 
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where C|TT>.I+) * s tne projection of the two-spin state Iff) 
onto the renormalized cell state |+), etc. The projection 
operation is a mapping from an Hilbert space of N spins, 
with dimension 2 N , to an Hilbert space of N/2 block 
spins, with dimension 2 N ' 2 . Thus f x : <C N h-> <£ n / 2 . We 
also define the embedding operator 



f 2 (a,a) = y^Cn^n \n) (n'\ 



(20) 



which may be cast as a 2 N x 2 N < 2 matrix whose columns 
are linear combinations of the right ground states of Hq. 
We have 



To = 



c l+).|tt> C |-),|TT> \ 

C I+>,IU> C |-),|T1) 

C I+),I1T) c h),|lT> 

V C I+},IU> c \-)Ml) / 



®JV/2 



(21) 



where now c\+\ I TT> * s ^ ne component of the renormal- 
ized cell state |+) that is embedded in the two-spin state 
|TT). Thus T 2 : £ N ' 2 h-> C N . We demand that if we 
apply the embedding operator followed by the projection 
operator we recover the identity on the 7V/2-dimensional 
block-spin space: T\T 2 — 1jv/2- But since the projection 
operation does not retain all the degrees of freedom of 
the system, projection followed by embedding does not 
yield the identity on real-spin space. Thus T 2 T\ ^ In- 
The renormalization prescription is then 



H'(a') =T l (a',<T)H(a)T 2 (a,a'), 



(22) 



where H' is a renormalized evolution operator. This pre- 
scription projects the original 2 N x 2^ Hamiltonian onto 
a 2 N / 2 x 2 N ' 2 subspace. If this subspace is suitably cho- 
sen, the renormalized Hamiltonian H' will have the same 
form as H, but with renormalized couplings A' = /(A) 
and rates V = p z T. From these relations one can deter- 
mine fixed points and critical exponents. 

We will now apply this scheme to the KCICs. 



right ground-state eigenvectors, (0,0, 0, 1) T and (c/(l — 
c), 1, 1, 0) T . We will use these to build T\ and T 2 , subject 
to the following constraints: 

1. The RG transformation must preserve probability 
conservation. Thus each column of (H FA )' must 
add up to zero. 

2. We require that V' a a+i has the same form as Ho, so 
that we can identify unambiguously the renormal- 
ized parameters. Note that by building T 2 from the 
ground states of Ho we ensure that the renormal- 
ized intra-cell Hamiltonian vanishes, i.e. H'o = 0. 

3. We must respect the fact that the FA model is triv- 
ially irreducible for all T ^ [l3j . This means that 
any configuration (bar that with all spins down) 
can be reached from an initial high-temperature 
configuration. This suggests that any two-spin 
state with at least one up-spin, namely |TT),|TI) 
and ITT): should be projected onto |+). 



4. Normalization. We require that T1T2 = 



One choice satisfying these criteria is 





( 



T 2 



2-c 



V 







2-c J 



(24) 

The matrix T\ projects Tt)JTI) an d lit) on t° |+) 5 
and |TT) onto |— ). T 2 embeds the state |+) as (2 — 
c)- 1 {c |TT> + (1 - c) |U) + (1 - c) |TT)}, and |-) as |||). 
The form of the ground state vectors for the FA model 
stipulates that the states |TT) an d |TT) are treated on 
equal footing during the projection and embedding op- 
erations, as befits a model whose dynamical rules are 
isotropic. We will see in the next section that this is not 
so for the East model. 

Using (J2U and we find 



V. RENORMALIZATION OF THE FA MODEL 

The FA model is defined by Equation (fl2*|l . When the 
blocking parameter p — 2, we need only consider the 
16 x 16 matrix 

H FA = (n<8*)<8(l®l) + (l®n)®(*®l) 

+ (£<g>n)<g> (1 ® 1) + (1 <g> <g> (n®l). (23) 

The brackets indicate the groupings of cells into blocks. 
The first and third terms in l)23|) comprise the intra-cell 
Hamiltonian H.q; the second and fourth terms are the 
inter-cell interaction V a>a +i- 

We must calculate the left and right ground states 
of the matrix Ho (|12fl . There are two left ground- 
state eigenvectors, (0,0,0,1) and (1,1,1,0), and two 



FA\i 





/ 


Z 2-c 


— c 


— c 


o\ 


1 




(1-cf 
2-c 


c 








2(2 - c) 




(1-cf 
2-c 





c 







V 











0/ 



(25) 



We can deduce the flow of the temparature parameter as 
follows. Let the ratio of the sum of the rates of the pro- 
cesses — > AA and A — > AA to the sum of the rates of 
AA — > %A and AA — > be A. Thus its unrenormalized 
or 'bare' value is 



An 



(1,2) + (1,3) 
(1,1) 



1 



= e- 1 ' T , 



(26) 



where (i,j) is element (i, j) of matrix (|12fl . Hence A — > 
as T — » 0. We can work out how A renormalizes by 



5 



7f 



a =pa 



C=- oo 



c=0 



c=l 



FIG. 1: FA model RG flow diagram for the temperature 
parameter c = (1 + e 1 ' T )~ 1 . The zero-temperature and high 
temperature fixed points, c* = (A* = 0) and (c* = 1) 
(Aa = oo), are respectively unstable and stable. The fixed 
point c* = -co (A* = —1) is unphysical, and inaccessible. 



calculating a similar ratio using the matrix The 
resulting RG recursion relation is 

A' = A(2 + A), (27) 

where A' is the renormalized counterpart of A. Equation 
(I27|l describes the flow of A away from an unstable zero- 
temperature critical point A* = 0, towards a stable high- 
temperature fixed point A* — > oo. The unphysical fixed 
point A* = — 1 is inaccessible. Figure ^ shows this flow. 

The RG procedure for the FA model using larger block 
sizes is unambiguous, because one obtains at each stage 
only two right- and two left ground state eigenvectors of 
the intra-block Hamiltonian. Thus for p = 3 we construct 
the following projection and embedding operators: 



T (3) = 



11111110 
1 



(28) 



and 



r 2 (3) =A/-(A) 



/A 2 
A 
A 
1 
A 
1 
1 



(29) 



V Af(A)- 1 J 

where Af(A) = (3 + 3A + A 2 ) -1 . In T 2 (3) , one inserts in 
the relevant slot one power of A for every up-spin in excess 
of one, in order to reflect the thermal suppression of these 
states. Thus the state | T T T ) (corresponding to element 
(1,1) of is penalized by a factor A 2 , wheras the state 
lit) (corresponding to element (5, 1) of (J2HJ) receives a 
penalty of one power of A. The generalization to larger 
block sizes is straightforward. We find that for general 
block size p the RG recursion relation is 



A fc = (l + A^xf-l, 



(30) 



where Afc is the value of A following the fc-th iteration 
of the RG. As expected (and required by the semi-group 
property of the renormalization group) we see from Equa- 
tions H27fl and (|30|l that two successive coarse-grainings 
using a block size of p = 2 is equivalent to one coarse- 
graining using a block size of p — 4. Thus A' = pA+C(A 2 ) 
near the critical point A* = 0. 



FIG. 2: An illustration of how length scales change under the 
RG blocking procedure. Say the system admits a dimension- 
ful correlation length £d . Then in the unrenormalized system 
(left panel) we may construct its dimensionless counterpart, 
measured in terms of the lattice spacing a: £ = £r>/ a - In the 
renormalized system (right panel) the new lattice parameter 
is a' = pa. Thus the renormalized dimensionless correlation 
length is £' = in/a' = 



The divergence of the dynamical correlation length fol- 
lows from standard RG arguments [23L |24| . Because the 
dimensionful correlation length must remain invariant 
under the RG transformation, the dimensionless corre- 
lation length £, measured in terms of the lattice spacing, 
must decrease by a factor of the blocking parameter, p: 
£' = p^ 1 ^ (see FigureEJ. We can write this relation as 



(31) 



In Equation (|31|l A' is the renormalized version of A. If 
we can write the RG equation for A near criticality in 
the form A' = p y X + C(A 2 ), then the correlation length 
is a function satisfying £,(p v X) — p^ 1 £,(X). Hence £(A) oc 
A -1 ^, where v± = 1/y. From Equation (|30|l we see that 
vj_ = 1, and hence near the critical point 



£(Ao) 



Ao 



(32) 



This correlation length corresponds to the characteris- 
tic spatial extent of structures ('bubbles') in space-time 
trajectories of the FA model at low temperature. We 
show one such trajectory in Figure [3] 

We can obtain the dynamical exponent z by noting 
that in the limit of zero temperature the nonvanishing 
elements of (|25|) are one-quarter those of (|12fl . We find 
that for general p the corresponding rescaling factor is 
p -2 . We interpret this factor as a rescaling of time under 
renormalization, defining the dynamical exponent z via 
t' = p~ z t. Thus for the FA model z = 2, signifying diffu- 
sive behaviour. This is as expected: the low-temperature 
dynamics of the FA model is known to proceed by diffu- 
sion of isolated defects 

We can infer the consequent relaxation time of the FA 
model by using the relationship between time and length 
scales, t ~ P, where I is the length scale being probed. 
Since the equilibrium length in the FA model scales as 
l eq ~ c -1 — see below, and Refs. — an d since the 

microscopic timescale goes as c, we expect the equilibra- 
tion time to have the leading order temperature depen- 
dence cr eq ~ c~ 2 T eq ~ exp(3/T). This scaling is 
known from previous work on the FA model [13J . 



One may also calculate 0] the density of excitated 
sites, n = |(1 + <r), both in the steady state and near 
the critical fixed point. The former is trivial for the 
FA model, since it obeys detailed balance, and one may 
therefore consider the calculation of the steady-state den- 
sity a test of the RG scheme. 

First note that the renormalization of the number op- 
erator does not depend on whether n sits in the left or 
right slot of the block: (1 ® n)' = (n <gi 1)' = (2 - c)~ 1 n a , 
for a block size p = 2. The RG recursion relation for the 
density then reads 



1 + Afc 

2 + X k 



(33) 



where the subscript k denotes the parameter obtained 
following k iterations of the RG. 

To extract the steady-state density we follow and 
write n(Afe) = a(Xk)n(Xk+i), where a(x) = (1 + x)/(2 + 
x). By iterating this equation along the RG flow we get 



n s (X ) 



\\»-X, 



i=0 



n(X* a ), 



(34) 



where n s is the steady-state density, and «(A*) is the 
density at the attractive fixed point A* = oo. Again fol- 
lowing we define G„(A) = YYi =0 a (Xi)- From Equa- 
tions l|27(l . (|33|l and the definition of a(X) we can write 



a (\ k \ - 1 d\nX k+ i 
2 d In Xk 

We can therefore write G„(A) as 

n ,n _ 1 d\nX n+1 
" 2 „+i dlnAo 



(35) 



(36) 



From l(2"T)) we have that A n +i = (1 + An) 2 ™ +1 — 1. Using 
this result with Equations (|34[) and (|36J1 . we get 



n s (Ao) = lim 



An 



oo Aq + 1 



A 



ra+1 



n(K)- (37) 



As n — » oo, A n +i — > oo, and so, noting that n(A*) = 1, 
we obtain the steady-state density 



n s {X ) 



A 



Ao + 1 



(38) 



This is as expected: detailed balance with respect to 
the Hamiltonian H(a) = ^J2j a j implies (rij) eg = 

E W [(! + <*) /2] e-W)J £ W = c. 

Near criticality we can write n(A) = p n(pA), and so 
n(An) ~ Aq. Thus the density vanishes close to criticality 
asn~c 8 with (3=1. 



100 



x 50 
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FIG. 3: Equilibrium space-time trajectories at T = 1.0 for 
the 1 + 1-dimensional FA (left) and East (right) models, re- 
produced from flU . Up-spins are black, down spins white. 
Space runs along the vertical direction, encompassing 10 5 
spins. Time runs along the horizontal. The characteristic 
length scales of both systems correspond to the vertical ex- 
tent of the 'bubbles' of down-spins, which scale in equilibrium 



The horizontal extent, r, of the bubbles is deter- 



mined by the relation r ~ l z , where the dynamical exponents 
for the FA and East models are z — 2 and z cx 1/T, re- 
spectively. As one observes the trajectories shown above on 
smaller length and shorter time scales, one moves from right 
to left along the RG flow diagrams shown in Figures 0] and Q 
Hence one eventually probes behaviour controlled by the crit- 
ical fixed point. For the East model, the disappearance of the 
characteristic length is consistent with the emergence of the 
fractal structure of the bubble boundaries. 



VI. RENORMALIZATION OF THE EAST 
MODEL 

The East model is defined by Equation . To renor- 
malize it using a blocking parameter p = 2, for example, 
we need only consider the 16 x 16 matrix 

H E = (n®£) ® (1 ® 1) + (1 <g> n) <g {i® 1). (39) 

The brackets indicate the groupings of cells into blocks. 
The first term in 139fl is the intra-cell component Ho; the 
second is the inter-cell interaction V a , a +i- 

We must calculate the left and right ground states of 
the matrix Ho dJ. The left ground states are rep- 
resented by the row vectors (0, 0, 0, 1), (0, 0, 1, 0) and 
(1,1,0,0). The right ground states correspond to 
the column vectors (0, 0, 0, 1) T , (0, 0, 1, 0) T and [c/(l - 
c), 1,0, 0] T . Next, we choose the projection and embed- 
ding matrices, which we call R\ and R2 so as not to con- 
fuse them with their FA model counterparts. One choice 
satisfying criteria 1-4 (see above) is 



Ri = 



1110 
1 



R 2 = 



I ac 
a(l-c) 
1 -a 

V 






1/ 



(40) 



where < a < 1 parameterizes a degree of freedom. This 
arises because the East model admits one more ground 
state vector than the FA model. Ri projects |TT) , |T|) 
and ||t) on t° an d III) onto |— ). i? 2 embeds the 
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state |+) as ac ] TT> + a(l - c) \U) + (1 - a) |4t>, and |-) 
as HI). With this choice we get 



{H E )' = {l + ac-a) 



( a(l-c) 2 -c 0\ 

-a(l - c) 2 c 



V 0/ 



(41) 



We deduce the flow of the temperature parameter in a 
similar way to before: let the ratio of the rates of the 
processes A% — > AA and AA — > A% be A. Then 



A 



(1,2) 
(1,1) 



-l/T 



1 



(42) 



where {i,j) is element of matrix (|llfl . Hence the 

bare temperature parameter has the same interpretation 
as in the FA model. We can work out how A renormalizcs 
by calculating the ratio of elements (1,2) and (1,1) of 
matrix (|41|l . The resulting RG recursion relation is 



A' = a _1 A(l + A), 



(43) 



implying an unstable zero-temperature critical point, 
A* = 0, as expected. 

The dynamical exponent z follows immediately. In the 
critical limit A — > 0, element (1, 1) of matrix 1)11(1 becomes 
unity. Hence we may interpret the renormalized value of 
this element as the time rescaling factor 2 _z . From (f 4 If) 
we get 



lim |l + a(c — l)|a, 



(44) 



and so z depends on the value we choose for a. 

Let us choose a. This parameter measures the extent to 
which we treat the states |lf) and on equal footing. 
In a model with symmetric dynamical rules, such as the 
FA model, we must treat these states identically. But the 
East model has asymmetric dynamical rules, suggesting 
that at some point in our calculation we must suppress 
||t) relative to |||), or vice- versa. At which point should 
we do this? We note that the projection matrix R\ treats 
and identically. If this were not the case, and 
we instead (for example) used 



R[ = 



110 
11 



(45) 



much less important to the dynamics than a configura- 
tion |TI), which is mobile. Thus when one renormal- 
izes the lattice using Ri and R2 with a = 1, one effec- 
tively discards dynamical pathways mediated by blocks of 
'jammed' spins |J.|). The RG process discards inaccessi- 
ble pathways in trajectory space {o~t t , o~t 2 ,•■■}, according 
to rules imposed by the Liouvillian of the dynamical pro- 
cess. Loosely, the projection matrix R\ identifies those 
single-spin states which are facilitating, whereas R2 picks 
out those two-spin states which are (internally) mobile. 

Setting a = 1 immediately yields a temperature depen- 
dent dynamic exponent: from Equation (|44|l we obtain 
2~ z = c, or z = (Tln2) _1 . Were a < 1, z would be 
independent of temperature to leading order. We thus 
conclude that maximal spatial anisotropy in the embed- 
ding process is a necessary condition for a temperature- 
dependent dynamic exponent. 

The RG scheme for the East model can be general- 
ized to larger block sizes. However, this procedure is less 
straightforward than for the FA model, because of the 
freedom one is afforded by the East model's many ground 
state eigenvectors. Furthermore, the results one obtains 
depends on whether one coarse-grains using a blocking 
parameter p equal to a power of 2, or not. 

Let us first illustrate the generalization of the proce- 
dure for the case p = 4 = 2 2 . We show that the results 
are consistent with the p — 2 scheme. We then argue 
that one should obtain a different dynamical exponent if 
one coarse-grains the system using a block size not equal 
to a power of 2, and then show explicitly for p = 3 that 
this is indeed the case. 

(v) 

Consider p — 2™, where n is an integer. Building R\ 
is straightforward: it is identical to t[ p \ its FA model 
counterpart. Thus R± is a 2 x 16 matrix whose top 
row is composed of Is apart from the rightmost element 
which is zero. Vice versa for the bottom row. 

The form of the embedding matrix is less obvious, be- 
cause the number of ground states increases as one in- 
creases the block size. However, we are guided by the 
form of the Liouvillian, which for block size p = 4 may 
be written schematically as 

TL E = (n<g>f<g>l(g)l) + (l<g>l<g)n<g>£) 

+ (1 <g> n <g> £ <g> 1) (46) 

+ (£<g> 1(8)1® 1) + (1 Ol® l®n). 



we would violate criteria 2 and 3 above. (R[ imposes a 
symmetry between flipping spins f <->J. in a two-spin block 
and flipping the resulting renormalized spin |+) |— ).) 
Therefore, we conclude that the embedding matrix Ti 
must treat and \\\) asymmetrically. The simplest 
way of doing this is to set a = 1, thus suppressing com- 
pletely the state If). This corresponds to the asser- 
tion that a spin configuration |||) (which is unable to 
change state unless connected to neighbouring spins) is 



Brackets again denote the grouping of cells into blocks. 
We take the first line of Equation l|46|l as the intra-cell 
Hamiltonian Tig. The second line vanishes under renor- 
malization as a consequence of T\ acting on it from the 
left, and so we ignore it; the third line comprises the inter- 
ccll interaction whose renormalization properties we wish 
to study. 

Guided by the form of R\ , we find that one choice of 
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satisfying criteria 1-4 above is 

/ A 3 
A 2 

A 
A 2 
A 

1 










i4 4) =A^(A) 



(47) 



V A^W(A)- 1 / 

where A/" (4) = (l + 2A + 2A 2 + A 3 )~\ We see that Equa- 
tion (|47|1 can be obtained from its FA model counterpart 
by using a simple rule-of-thumb: suppress all states of 
the form || • • • ) [corresponding to elements (9, 1)— (18, 1) 
of (|4*7Jl ] , as well as states possessing a 'frozen' up-spin at 
the right-hand boundary of the block. Thus states |TT4T) 
[corresponding to element (3, 1) of l|4*7jl ] and [Tilt) [ele- 
ment (7, 1)] have been removed. We see again that the 
embedding operator plays the role of a dynamical 'filter', 
eliminating those states which play a sub-dominant role 
in the dynamics of the East model. 

With these choices of embedding and projection oper- 
ators we obtain the RG recursion relation for the tem- 
perature, 



A' (4) =A(1 + 2A + 2A 2 + A 3 ), 
and a relation for the dynamical exponent: 



= lim 



A 2 (l + A) 



a^o (l + A + A 2 )' 



(48) 



(49) 



Equation (|48() is identical to the result one would ob- 
tain via two coarse-grainings using a block size p = 2 
[Equation (|43|l ]. as required. Equation l|49|l yields the 
dynamical exponent z — (Tin 2) , as before. 

We shall demonstrate how one can generalize this ap- 
proach to arbitrarily large n. Let us use reaction-diffusion 
notation (f— » 1, 0), and write the projection and em- 
bedding operators in the form 

4 = |+)(e<i*M +|->(oo---o|, (50) 



and 



B.2 = ^a u \l*)J 



|00 • • ■ 0) 



(51) 



The symbol 1* denotes any state 1 1 ■ ■ ■ ) starting with a 
1, and the {a^} are a set of coefficients [see e.g. the first 
column of l|47Jl]: the projection operator R\ allows one 
to compare this notation to the matrix representations 
employed previously. The normalization requirement R\ ■ 
R.2 = 1 implies J2i* a i* = 1- Thus at least one of the 
coefficients a must be of 0(1). 

The values of these coefficients are fixed by the eigen- 
vectors of the intra-block evolution operator, as we have 
discussed. We can see how these coefficients deter- 
mine the properties of the model under renormaliza- 
tion, as follows. We find that 'bulk' states of the form 
(•••l(E>n(E>i®l---) vanish under renormalization as 
a consequence of the projection operator acting from 
the left. We are therefore left with the 'surface' terms 
(£<g> 1 (g) • • • <g> 1) and (1 Cg) • • • ® n), in which the operators 
n and I sit at the edge of the block. We find that, under 
renormalization, 



Ri ■ (l 



h)-R 2 




(52) 



and 



Ri 



l\-Ri 



(oio-o) % 



A 



Oio-.-O 



(53) 

In Equations (|52() and l|53|) primes denote renormalized 
operators. The symbol 1*1 denotes states |1 • • • 1) start- 
ing and ending with a 1 . In (|53|) the temperature param- 
eter A has been rescaled by the coefficient aio-.-o which 
weights thermally the state 1 10 • • • 0) with a single 1 at 
the leftmost edge, followed by a string of 0s. From the 
previous discussion we know that this coefficient is of or- 
der unity, and hence the recursion relation for A will be 
marginal, as we have found. 

The dynamical exponent follows by noting that the 
product of Equations l|52|l and (|53|l constitutes the renor- 
malized evolution operator, and so the prefactor de- 
scribes the rescaling of time as a consequence of rescaling 
space. Thus 



oc lim (etin. 



(54) 



The constant of proportionality in l|54|l is (1 — c) -1 = 
1 + A, i.e. the reciprocal of element (1,1) in the unrenor- 
malized East model Hamiltonian l|ll|l . The first factor 
on the right hand side of l|54|> is of order unity. The 
second factor is fixed by the embedding operator, which 
is in turn determined by the relevant East model eigen- 
vectors. The rule-of-thumb we obtained above tells us 
that we remove from this factor any state with a frozen 
rightmost up-spin. This may be regarded as an entropic 
suppression of states playing only a sub-dominant role in 
the dynamics. Those states starting and ending with a 1 
which are important for the dynamics of the East model 
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are for block sizes p — 2, 4, 8 and 16, 

11 (A) , 1011 (A 2 ) , 10001011 (A 3 ) , 

1000000010001011 (A 4 ) . (55) 

All have a 'mobile' rightmost up-spin. The thermal 
weighting of each state is given in brackets. These states 
are important because of the hierarchical dynamics of the 
East model 0, El i which dictates that two defects sep- 
arated by a distance d are relaxed by establishing a set 
of isolated defects between them, at distances d/2, 3<i/4 
etc. Thus for block size p = 4 the dominant dynami- 
cal pathway proceeds via the state 1011, with a thermal 
weighting of A 2 (and not, for example, 1001, which has 
a weighting of A). Hence axon ~ A 2 , aiooi = 0, and 
so lim>^o a i*i ~ Consequently, the rescal- 

ing factor is A~ z ~ A 2 , and the dynamical exponent 
z = (Tln2) _1 , as required. (In the FA model, states 
such as 0001) are permitted, leading to temperature- 
independent a coefficients and hence to a temperature- 
independent z.) Thus i?2, which attempts to reconstitute 
an unrenormalized state from a coarse-grained state, cap- 
tures both energetic effects (the powers of A weighting 
thermally the various states) and entropic effects (the 
'zero' entries corresponding to those suppressed entrop- 
ically). We conclude that the RG scheme for the East 
model generalizes readily to larger block sizes. 

It is interesting to note that if one uses blocks of size 
not equal to a power of 2, one obtains a slightly different 
result for the dynamical exponent. We argue that this is 
a consequence of the hierarchical dynamics of the East 
model takingplace naturally in blocks of lengths equal to 
a power of 2 ■ We can derive the approximate value 

of z that one should obtain from a coarse-graining over 
block sizes p ^ 2". Let us take p = 3 as an illustration. 

Consider the coarse-grained relaxation process |++) -^-» 
H — ). We wish to determine the leading order tempera- 
ture dependence of 7', noting that the rate for the equiva- 
lent unrenormalized process, Iff) |ff), is 70 = 0(1). 
From our previous discussion of the form of the embed- 
ding matrices, we can infer that the dominant dynamical 
pathway (involving 'unrenormalized' spins) contributing 

to this renormalized process is Iff If If) Iff HID- 
To relax the second up-spin, one must create two ex- 
tra up-spins to the right of the first up-spin. Hence this 
pathway has a rate 71 ~ c 2 , and the renormalized rate 
0(c 2 ). Other pathways also contribute to the renor- 



1 — ), but do so either with 
If If III)) — in which case 7' 



r 

malized process |H — (-) — 
rates ~ c 2 (e.g. Iffffff) 
is changed only by a temperature-independent numerical 
factor — or with rates higher order in c (e.g. the pathway 
Iffffff) — > Iffffff))- These we may ignore. Since we 
interpret the overall rescaling of the fundamental relax- 
ation rate deriving from a coarse-graining of space as the 
numerical factor p~ z , we would therefore expect for p = 3 
that 3~ z ~ c 2 , or z ~ 2/(Tln3). 

Loosely, then, we expect that by coarse-graining space 
in blocks of size 2"~ 1 < p < 2" one should obtain z w 



n/ (Tlnp) [which tends to z — > (Tln2) _1 when p — > 
00] . Coarse-graining using block sizes p — 2" yields z — 
(T In 2) _1 . This is as we expect: the energetic barriers for 
relaxing chains of lengths 2" _1 < p < 2™ and p — 2™ are 
identical, but the entropic barriers are larger for the latter 
case. Thus one would expect the dynamical exponents 
to differ. More sophisticated arguments [2!| reveal that 
z is bounded by (Tln2)- 1 and (2Tln2)" 1 . 

We can show that our guess for the dynamical expo- 
nent is borne out in the case p — 3 by the RG scheme. 
We construct the embedding operator according to the 



i?f=Af (3) (A) 



/ A 2 
A 

1 






(56) 



V A^(A)- 1 J 

where Af^ = (l + A + A 2 ) 1 . Together with the obvious 
choice for the projection operator we find 



lim 

(l + A + A 2 )' 



(57) 



yielding z = 2/ (Tin 3), as advertised. We conclude that 
the RG scheme for the East model can be generalized 
to larger block sizes, but more naturally so for the case 
of a blocking parameter p equal to a power of 2. For 
simplicity we shall focus on the case p — 2. 

With a = 1 the RG recursion relation (|43|l may be 
iterated to give 



A p = A + n\ 2 Q + O (A 3 ) 



(58) 



where \ p is the value of the temperature parameter fol- 
lowing a coarse-graining of the system by a factor p = 2™. 
Since the bare value of A > 0, we see that de- 
scribes a system with an unstable zero-temperature crit- 
ical point A* = 0, and a stable high-temperature fixed 
point A* — > 00. Now, however, the temperature parame- 
ter A is marginally relevant near the fixed point A* = 0. 
The RG flow diagram is shown in Figure 

To determine the correlation length in the East model 
we proceed as follows. From the recursion relation l|58|l 
we see that the correlation length satisfies 



x+ — x J = p 



(59) 



For small values of x = A, corresponding to low tempera- 
tures, we have no solution £(A) of (|5^|l to first order in A. 
Thus near criticality the East model possesses no char- 
acteristic length scale. This is consistent with the nature 
of the space-time trajectories seen in numerical simula- 
tions, such as that shown in the right panel of Figure 
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c=0 



c=l 



FIG. 4: East model RG flow diagram for the temperature 
paramter c = (1 + e 1 ^ T )~ 1 . The zero-temperature critical 
fixed point c* = (A* = 0) is unstable. Near this fixed 
point, corresponding to observations of the system on small 
length and time scales, the East model admits no charac- 
teristic length. This is consistent with the fractal structure 
of space-time trajectories seen in numerical simulations, such 
as those shown in Figure [3] The dynamical exponent z is 
proportional to 1/T, indicating a rapid slowing-down of the 
dynamics near the critical point. As one observes the sys- 
tem on progressively larger length and time scales, one sees 
the emergence of a characteristic length growing as the dou- 
ble exponential of reciprocal temperature. This length sat- 
urates rapidly, and is eventually supplanted as the charac- 
teristic length of the system by by the equilibrium domain 
size, whose properties are controlled by the high-temperature 
stable fixed point c* = 1 (A* = oo). 



These display a fractal structure ^jj , and hence possess 
no characteristic length. 

On sufficiently large length and time scales the sys- 
tem will reach equilibrium, at which point the heights of 
bubbles will be determined by the equilibrium spin dis- 
tribution. This does have a characteristic length. One 
therefore expects to see, sufficiently far from criticality, 
the emergence of a length scale. Below, we show that 
Equation (|59|l indeed admits a growing length in such a 
regime. This corresponds to the eventual 'blurring out' of 
the fractal boundaries of clusters as one observes the sys- 
tem on progressively larger scales. The emerging length 
scale corresponds to the spatial extent of bubble regions. 

We can quantify the emergence of this length by con- 
sidering an infinitesimal RG transformation. The block- 
ing parameter p is necessarily an integer, because our 
model is defined on a lattice. But we can generalize p by 
considering an infinitesimal change of scale according to 
p = 1 +1, where I <C 1. By writing £(A + £A 2 /ln2) « 
£(A) + £d£(\) /dl and X' - A « IdX/dl we obtain the flow 
equations for the temperature and correlation length: 



and 



d\{£) 
dl 



dm 

dl 



X{1) 2 / In 2, 



= -m- 



(60) 



(61) 



The intial data for Equations Ij60|l and (|61[) are X(£q) = Ao 
and £(lo) = £o ; respectively, where the subscript zero de- 
notes an unrcnormalizcd (physically meaningful) quan- 
tity. The parameter £q acts as a short distance regulator 
(or ultraviolet cutoff), and should be taken to zero at the 
end of the calculation. 

One now iterates the RG by integrating l|60|l until 
A(f) = 0(1), yielding I - l « In2/A . From we 



obtain £(£) = £oe ^ l °\ and so the correlation length 
varies with temperature according to 

Co(A ) ~ exp (1/(A ln2)) ~ exp (e 1 ^/ In 2) . (62) 

Away from the critical point A* = we therefore see 
an extremely rapid growth of the dynamical length scale 
with temperature. 

This length scale corresponds to the emergence of a 
characteristic length ^ away from criticality, and not 
to an equilibrium length scale l eq . The latter may be 
defined as the reciprocal of the particle density in the 
steady state (see below), and scales as c~ x . The dynam- 
ical length is a nonequilibrium critical quantity, and will 
be cut off rapidly as one probes larger length and time 
scales. Thus, in terms of the RG flow diagram, Figure 01 
the steady-state behaviour is obtained near the attrac- 
tive fixed point A* — > 00, where one probes length and 
time scales much larger than those on which critical fluc- 
tuations are manifest. The critical behaviour will be ob- 
served on short length and time scales, near the critical 
fixed point A* = 0. 

The characteristic equilibration time follows from the 
relation r ~ l z , where I is a typical length scale. We 
have z = (TTn2) _1 . Taking the equilibrium domain 
length l eq ~ A -1 , we find the equilibration time scale 
Teq ~ c ~ z ~ A lnA = exp{l/(T 2 ln2]|. This agrees with 
results obtained by other means [Tj. We assume that 
the dynamical exponent z obtained near criticality holds 
in the region of the attractive fixed point. 

One may also calculate the density of particles, 
n = ^(l + <j). First note that the number operator 
renormalizes differently depending on whether n sits in 
the left or right slot of the block: (n (g) 1)' = n a , versus 
(lCgm)' = cn a . Hence we will define our density operator 
as n — \ (n ® 1 + 1 (8> n). The RG recursion relation for 
the density then reads 



1 / 1 + 2A* 



n k = - 



1 + A A 



(63) 



To extract the steady-state density we write (|63|l as 
n(Afc) = a(\ k )n(\ k+ i), where a(x) = (l + 2x)/ [2(1 + x)]. 
By iterating this equation along the RG flow we get 



n s {X ) 



i=0 



(64) 



where n s is the steady-state density, and n(A*) is the 
density at the attractive fixed point A* = 00. Next, define 
G„(A) = IlLo a ( A «)- From Equations ® and (ESJ we 
can write 



ld\n\ k+1 

s At = 

1 ' 2 dlnAt 



as in the FA model. Hence 



G, 



A,, 



(x) = T ^-(i + o(x-l 1 )). 

Ao + 1 



(65) 



(66) 
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Taking n — > oo gives A„+i — > oo, and by noting that 
n(A*) = 1 we obtain the steady-state density 



n s (X ) 



Ao + 1 



(67) 



as expected. 

The behaviour of the density near criticality (A* = 0) 
follows from the relation n(AA) = a(X*)n(X'), where 
AA = A — A* = A. If we iterate the RG until the renor- 
malized coupling A' = 0(1), i.e. £/£ ~ e ln2 / A °, wc find 
n(Ao) ~ e~ ln2 Ao . Thus the density vanishes close to the 
critical point faster than any power of T. 



VII. RENORMALIZATION OF THE BCIC 

In this section we will apply the RG scheme to the 
BCIC, a model whose kinetic constraint interpolates be- 
tween that of the East and FA models. We find that on 
suitably large length and time scales (or for suitably low 
temperatures) the BCIC behaves like the FA model. This 
agrees with existing numerical and analytical results |20| . 

The ground state eigenvectors of the BCIC (llUf) are 
the same as those of the FA model. If we use jSU and 
|[TDJ we find 



/ 2=S - c -be -be \ 



(H 



by 



V 



b 
b 

(c-2) " 





be 
h be 



be 
be 







0/ 



(68) 



Equations (|10fl and 168|) yield the same recursion relation 
for the temperature parameter A as in the FA model, 
Afe+i = Aa;(2 + Afc). They also yield a recursion relation 
for the asymmetery parameter b: b k +i = b k . Thus the 
asymmetry b is a marginal operator, and does not flow 
under renormalization. From the RG relation for A, we 
see that for any b £ (0, 1) the interpolation model falls in 
the universality class of the FA model, rather than the 
East model. 

However, we expect the interpolation model for small 
values of b to display a crossover from East- like to FA-like 
behaviour [2(J. This suggests that by projecting H b onto 
a subspace spanned by only the ground states of (fTU|) we 
have omitted this crossover behaviour. We can recover it 
in the following way. 

First, we note that the difference between the East and 
FA models manifests itself in the treatment of the states 
|T J.) and |TT) during embedding. In the East model the 
latter is completely suppressed (see Equation 1(10)1): in the 
FA model, both are treated on equal footing [Eq. 1|24|) ]. 
By restricting our RG scheme to a subspace of the ground 
states of Equation ijTUI) . we are unable to construct an 
embedding operator that treats ||T) an d IT!) asymmet- 
rically [cf. i? 2 , Equation (|3U|) ]. 

To remedy this, we now include the first excited right 
eigenvector of (|l(Jfl in our embedding operator. We will 



call this eigenvector e. This is akin to calculating higher- 
order "loop" diagrams to check if b, ostensibly a marginal 
operator, is relevant at second order, e has eigenvalue 
2(1 — b)bc + 0(c 2 ), and is therefore a "gapless exci- 
tation" in the East model limit, 6^0. Note that 
e = (ei,e2, 1,0) T , where the are functions of c and b. 
For small c we have e « |(2& - l)c, -1 + (1 - 2b)c, 1, 0]. 
Let us now construct a new embedding operator, 

/ c + ae\ \ 
1 — c + ae2 
1-c+a 

2-c 



R* = 



(69) 



and demand that in the limits b 



ts o — > and b — > | we 
recover the respective embedding operators for the East 
and FA models, namely i?2 and T^. This is achieved by 
setting a = -(1 - 26) (1 - c). We note that R X R 2 = 1. 

Our renormalization prescription is now (7i b )' = 
R\H b R2- We derive recursion relations for the param- 
eters c and b in a similar way to before: we define the 
unrenormalized temperature parameter A as the ratio 



(1,2) + (1,3) 



1 



(70) 



where (i,j) is element of the matrix Tib = C, Equa- 
tion (|10|) . We define the renormalized parameter A' by 
the ratio of the corresponding elements of the renormal- 
ized matrix (7i b )'. This gives us the recursion relation 
Afe+i = /(Afc,/ife). The parameter /i is the scaled asym- 
metry parameter, whose unrenormalized value we define 
as 



/' 



(1,3) b 
(1,2) 1-6" 



(71) 



The elements (i, j) again refer to Equation l|l(J|l . We write 
the recursion relation for //, obtained from the elements 
of (H b Y, as fik+i = 9(Afe,/ifc). 

The behaviour of the functions / and g thus determine 
the crossover properties of our model. We find that A has 
an unstable zero-temperature fixed point A* = 0, and an 
attractive high-temperature fixed point A* — » oo. The 
asymmetry fj, has an unstable maximal-aysymmetry fixed 
point /x* = 0, corresponding to the East model, and an 
attractive symmetric fixed point fx* = 1, corresponding 
to the FA model. Thus any BCIC with less than maximal 
asymmetry will behave at long length and time scales like 
the FA model. FigureElshows the qualitative RG flow of 
the BCIC. 

For the case of p — 2 we find that 



A*(l + Afc)+/i(A fe )/i fe + 0(/4), Mfc~0; 



k+1 ~ \ A fe (2 + A fe ) + h{\ k )ixl + O (Al) , fJLk « 1, 

(72) 

where fi{x) = (2 + 9x + llx 2 + 6x 3 + x 4 )(2 + 3x + x 2 )- 1 , 
1)(2 + x)/[8(l + x)], and p, k = 1 - Li k - 



x 2 (x 
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Equations l|72[l thus reproduce the recursion relations for 
A in the East model and FA model limits, respectively 
Equations (|43|1 and (|27|l . The asymmetry parameter fi is 
a relevant perturbation, whose flow is governed by 



(i + ^) M)b -/3(Afc)/4 + e>(/4), 
■iik + o (Afe) , 



l+A, 



0; 
1, 

(73) 

where f 3 (x) = (2 + Ax + 5x 2 + 3x 3 )/[x 2 (l + x){2 + x)]. 

We can deduce the flow of the BCIC away from max- 
imal asymmetry, /i = 0, by studying Equations l|72l) and 
|JZ3J in the regime (j, < A < 1. Writing /3 A = A' - A w 
(AX/ di, and a similar relation for /x, we obtain 



£A^ 
A " 



A 2 



A1 + 0(A M )+0( M 2 /A); (74) 
0( M 3 /A 3 ) + 0(M 2 /A). (75) 



Equations (|74|) and l|75|) may be solved in terms of the ex- 
ponential integral function Ei (A -1 ), although the physi- 
cal interpretation of this solution is not obvious. We can 
more clearly determine the essence of the crossover as 
follows. 

The temperature parameter A has RG eigenvalue 
(East model) or 1 (FA model). It therefore grows much 
less rapidly than the asymmetry parameter fi, which has 
(initial) eigenvalue Aq" 1 ^> 1. Hence from (17311 we have 
jj = 2 y ^/i w Xq 1 ^, giving the RG eigenvalue for the 
asymmetry parameter as « (Tin 2) . Let us now 
write a standard RG scaling form for the particle den- 
sity, 

n fl (AV) =pn(py>\,py^,p- 1 tp- z t) (76) 

To derive a crossover temperature, we iterate the RG un- 
til p Vx = 0(1). The fi— dependent scaling combination 
is then A _ ^' J ' A . When this becomes large, i.e. 0(1), 
one would expect the BCIC to behave like the FA model. 
Taking for simplicity y\ = 1, we find a crossover temper- 
ature T xo ~ (— ln/i) . This scaling agrees with that 
obtained by equating the relaxation timescale for the [i- 
suppressed symmetric process, t$ ~ (/xA) with that for 
the asymmetric process, ta ~ exp (l/T 2 ln2) [TTL l2p| . 

We can extract crossover time- and length-scales from 
Equation l|76l) by iterating the RG until, respectively, 
= 0(1) and p~ z t = 0(1). These give £ xo - ^~ Tln2 



The real-space RG therefore confirms that for any- 
thing less than maximal asymmetry, the BCIC will on 
long length and time scales display FA-like, as opposed 
to East-like behaviour |2CI. 



VIII. CONCLUSIONS 

We have used the simple real-space RG scheme of 
Refs. 0S2 to derive the zero-temperature critical be- 
haviour of the FA, East and BCIC models. Our findings 




and t n 



-2Tln2 



FIG. 5: RG flow diagram for the BCIC, in (A, /x) space, /i = 
(resp. 1) corresponds to East (resp. FA) model behaviour. 
The critical fixed point (0, 0) is unstable; the attractive fixed 
point (oo, 1) corresponds to the high-temperature fixed point 
of the FA model. 



agree with known results 0, U El El H3 > but offer a 
different and unified approach to these systems. We are 
also aware of alternative real-space RG studies of KCICs 

mill. 

The real-space RG scheme used in this paper is suffi- 
ciently flexible to be extended to more complicated mod- 
els. An interesting possibility would be to use this scheme 
to study a recently-introduced model of the reentrant 
glass transition in colloids |28| . which combines dynami- 
cal constraints with static interactions. 
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